Synthetic spin-orbit coupling in ultracold A-type atoms 
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We consider the simulation of non-abelian gauge potentials in ultracold atom systems with atom-field inter- 
action in the A configuration where two internal states of an atom are coupled to a third common one with a 
detuning. We find the simulated non-abelian gauge potentials can have the same structures as those simulated 
in the tripod configuration if we parameterize Rabi frequencies properly, which means we can design spin-orbit 
coupling simulation schemes based on those proposed in the tripod configuration. We show the simulated spin- 
orbit coupling in the A configuration can only be of a form similar to p x o y even when the Rabi frequencies are 
not much smaller than the detuning. 
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I. INTRODUCTION 

Many interesting quantum phenomena have been found in 
condensed matter physics when electrons have a spin-orbit 
coupling, such as the spin Hall effect and the topological in- 
sulator 111, IH. Ultracold atom systems are now regarded as 
simulation platforms to study condensed matter physics [3], 
therefore it is important to realize spin-orbit couplings in these 
systems. It has been shown theoretically that abelian and non- 
abelian gauge potentials can be simulated for ultracold atoms 
via their interaction with laser fields, and non-abelian gauge 
potentials can be used to generate spin-orbit couplings [4- 
6fl. A general spin-orbit coupling for ultracold atoms has not 
yet been realized experimentally, however some experimental 
progress towards this direction have been made 071 41311 . 

Theoretical schemes to simulate spin-orbit couplings for ul- 
tracold atoms are usually proposed with atom-field interaction 
in the so called tripod configuration, where two dark states are 
used to form the effective spin space II 14142011 . These schemes 
have a drawback that the two dark states are not the lowest 
energy dressed states, hence atom-atom interactions can in- 
duce collisional decay. Some authors also consider the sim- 
ulation of spin-orbit coupling with atom-field interaction in 
the A configuration, where two lowest energy dressed states 
are used to form the effective spin space ll2ll I22I1 . Recently 
several experiments have realized the special spin-orbit cou- 
pling p x a y for ultracold atoms via Raman process lfTll - tl3ll . 
which is a scheme of A configuration with Rabi frequencies 
much smaller than the detuning. Some other kinds of methods 
are also proposed to simulate spin-orbit couplings in ultracold 
atom systems [23142511 . 

Although many interesting features of ultracold atoms have 
been theoretically found when they have a general spin-orbit 
coupling such as the Rashba and Dresselhaus spin-orbit cou- 
plings, currently we can only experimentally achieve the spe- 
cial spin-orbit coupling p x a y for ultracold atoms via Raman 
process. Since Raman process is a scheme of A configuration 
with Rabi frequencies much smaller than the detuning, it is 
natural to ask whether more general spin-orbit couplings can 



be simulated in A configuration when Rabi frequencies are not 
much smaller than the detuning. We find the answer is NO at 
least in our concerned A configuration. 

The structure of the paper is as follows. We first give an 
analytical expression of the simulated non-abelian gauge po- 
tentials in our concerned A configuration, which can help us to 
design spin-orbit coupling simulation schemes based on those 
proposed in the tripod configuration. We then consider a spin- 
orbit coupling simulation scheme where two plane waves are 
used in the A configuration. We find the simulated spin-orbit 
coupling can only be of a form similar to p x a y due to the non- 
degeneracy of the two lowest energy dressed states. We also 
analyze how the relative magnitude of the two lasers affect the 
simulated spin-orbit coupling Hamiltonian in this scheme. 



II. NON-ABELIAN GAUGE POTENTIAL SIMULATION IN 
A CONFIGURATION 



A general theory on the simulation of non-abelian gauge 
potentials for ultracold atoms is presented in Ref . J5J . Here 
we focus on atom-field interaction in the A configuration. As 
shown in FIG. 1, suppose two internal states |1) and |2) of an 
atom are coupled to a third common one |3) via laser fields 
with a detuning. The atom Hamiltonian will be 



H 



(1) 



where Hq is the atom-field interaction and V is the possible 
external trapping potential. In the interaction picture, 



H = ftA|3)(3| +ft[fii|3)(l| + fi 2 |3)(2| + H., 



(2) 
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where A is the detuning that is assumed to be positive, fij and 
0,2 are Rabi frequencies. We parameterize two Rabi frequen- 
cies as 



Cli = ^ tan 26» cos 0e iSl , fi 2 = — tan 29 sin 4>e lS2 , (3) 
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FIG. 1 : (a) The tripod configuration and (b) A configuration. If in 
the tripod configuration Rabi frequencies are parameterized as fli = 
f2 sin0 cos cj>e iSl , Q 2 = H sin 6 sin 4>e iS2 , and H 3 = A cos 6>e iS3 , 
and in the A configuration Rabi frequencies are parameterized as 
0,1 = ftan26>cos< ? !>e i(Sl - s ' 3) and Q 2 = § tan 26 sin tf>e i{S2 - S3) 
with A being the detuning, then the simulated non-abelian gauge po- 
tentials in the A configuration have exactly the same structures as 
those simulated in the tripod configuration. 



with — 7r/4 < 9 < tt/4. Under this parameterization the three 
eigenstates of Hq are 

|ei) =sin0e" aSl |l) -cos0e~ lS2 |2), (4) 
|e 2 ) =cos6»(cos0e- tSl |l)+sin< ? !)e~ 4S2 |2)) -sin0|3), 
|e 3 ) =sin6'(cos< ? !)e"' i ' Sl |l)+sin0e" lS2 |2)) +cos6»|3), 

with the coiTesponding eigenvalues being 



Ex = 0,E 2 = -HA 



cos(2(9) 



E* = hA 



cos(26>) ' 



(5) 



The full quantum state of the atom can be written in the 
form |$) = Yli=i *&i\ e i)- Due to the position dependence 
of the dressed states |e^), when we substitute the full quan- 
tum state |$) into the Schrodinger equation with H given in 
Eq. ([TJ, we can find that the column vector of wave functions 
^ = (vPi, ^>2jJi/3) T satisfies the Schrodinger equation with 
Hamiltonian ||5t] 



2m 



(6) 



where A and V are 3x3 matrices: 



A,, 



ih(e n \\Ve m ),V n 



E n 5 + {e n \V\e m ). (7) 



We note that A and V give contribution to the off-diagonal 
elements and E n give contribution to the diagonal elements 
of H. 

Since A > and -7r/4 < 9 < 7r/4, there are E 3 > 
Ei > E2 and E3 — E\ > HA. Now we discuss the con- 
dition when the off-diagonal elements of H connecting \I r 1 
to ^3 and ^2 to "J 3 can be neglected, so that it can be used 
to simulate the movement of a particle with spin- 1/2. Note 
that A n:Tn usually has a magnitude of momentum Pl of the 
applied laser fields, therefore the off-diagonal elements of H 
connecting ^1 to ^3 and \E , 2 to ^3 can be neglected when 



*C hA, |(e n |V|e m )| <C hA and atoms move very slowly 
(i.e., @- <C hA). If these conditions are satisfied the wave 
functions ( x l/i,Vl / 2) T will be approximately decoupled from 
^3 and evolute under the Hamiltonian ^ 



Q-A) 2 - 
2m 



(8) 



Here A, V and $ are 2 x 2 matrices, the elements of A and V 
are described in Eq. (0 with n,m = 1,2, and 



^n,m — A n 3 ■ A3 m , 

2m 



n, m = 1, 2. 



(9) 



The Hamiltonian H e f f in Eq. (O simulates the movement of 
a particle with spin- 1/2 in gauge potentials, where two low- 
est energy dressed states |ei) and |e2) represent spin up and 
spin down respectively. Here we emphasize that the 9 is not 
required to be small to get H e ff, i.e., the magnitudes of Rabi 
frequencies Oi and H2 are not required to be much smaller 
than the detuning A. 

We can get an analytical expression of the simulated gauge 
potentials A and $ in i? e // by substituting Eq. © into 
Eq. (0. But note that the dressed states |ei), |e 2 ) and ^3) 
have the same mathematical structures as the two dark states 
\Di), \D 2 ) H and the bright state \B) = \D ) [26] in the 
tripod configuration respectively, and in the tripod config- 
uration the simulated gauge potentials can be expressed as 
A n ,m = ih(D n \\VD m ) and = ^A nfi ■ A 0>m JH], we 
can conclude immediately that the simulated gauge potentials 
in our concerned A configuration have the same mathematical 
structures as those simulated in the tripod configuration, i.e., 
we can obtain an analytical expression for the simulated gauge 
potentials A and $ of H e f f in our concerned A configuration 
just through replacing S13 and S23 in Eq. (13) and Eq. (14) of 
Ref. |5[] by Si and S 2 respectively. Here we write them down 
for completeness : 

A u = h{cos 2 (j)VS2 + sin 2 <t>VSi), (10) 
A ia = Kcos6[-sin(2<i>)(?7 Si -VS 2 ) -iV<t>], 
A2.2 = hcos 2 9(cos 2 cj)\7Si + sin 2 <f>VS 2 ), 

$1,1 = ^szn 2 9[lsin 2 (2(b)(VS 1 - VS 2 ) 2 + (V0) 2 ], 
2m 4 

h 2 1 

$1.2 = — sm0[-am(20)(VSi-VS 2 )-iV<fl 
2m 2 

■[^sin(29){cos 2 (f)\7S 1 + sin 2 cjNS 2 ) - iV9], 
h 2 1 

$ 2 2 = — \-sin 2 (29)(cos 2 (j)VS 1 + sin 2 ' cjN S 2 ) 2 + (V6») 2 ]. 
2m 4 

We will give explicit examples of Si and S2 to show how 
spin-orbit coupling Hamiltonian can be obtained. 

Our parameterization of the two Rabi frequencies not only 
gives us a convenient way to obtain an analytical expression 
of the simulated gauge potentials as shown above, but also 
gives us a way to design spin-orbit simulation schemes in the 



3 



(a) 



|3> 



(b) 




FIG. 2: (a) Two plane waves coupling two lower internal states of an 
atom to a common higher internal states, (b) The directions of the 
two plane waves. 



A configuration based on those proposed in the tripod config- 
uration. Suppose there is a scheme to simulate spin-orbit cou- 
plings in the tripod configuration with the Rabi frequencies 
Cli = Ct sin 9 cos (f>e 1 , 2 = ^ sin 9 sin 4>e 2 and f2 3 = 
f2cos#e 3 , then in our concerned A configuration we can 
design a counterpart spin-orbit coupling simulation scheme 
using the Rabi frequencies Sli = — tan2#cos</>e 4 ( Sl ~ s ' 3 - ) 
and VL 2 = y tan 20 sin 4>e i( - S2 ~ S3 '> with A being the detun- 
ing. These two schemes simulate the same kind of spin-orbit 
coupling since the gauge potentials simulated by them have 
exactly the same mathematical structures, where the differ- 
ence is that there are some restrictions on the parameter 9 and 
there is a Zeeman term in H e ff due to the non-degeneracy of 
the two lowest energy dressed states in the A configuration. 
The known schemes to simulate spin-orbit couplings of Refs. 
I2H E2I1 in the A configuration can be regarded as the coun- 
terpart schemes of Refs. Il4l[l5ll in the tripod configuration 
where standing waves are used. In the next section we will 
consider a new spin-orbit coupling simulation scheme where 
only two plane waves are used, which can be regarded as the 
counterpart scheme of Refs. 116142011 in the tripod configura- 
tion where only plane waves are used. 



III. SPIN-ORBIT COUPLING SIMULATION IN A 
CONFIGURATION 

In this section we consider a spin-orbit coupling simulation 
scheme where two plane waves are used in the A configu- 
ration. This is an example to show the simulated spin-orbit 
coupling in the A configuration can only be of a form sim- 
ilar to p x a y even when the Rabi frequencies are not much 
smaller than the detuning. We also analyze how the relative 
magnitude of the two lasers affect the form of the simulated 
spin-orbit coupling Hamiltonian in this scheme. 

As shown in FIG. 2 two internal states |1) and |2) of an 
atom are coupled to a third common one |3) via two plane 
waves respectively with the same detuning. The directions of 
the two waves are in the x-z plane, and the phases of the two 
Rabi frequencies are assumed to be 

51 = ki ■ r = kxsmtp + kzcostp, (11) 

52 = &2 ■ T= — kxs'mip + kzcosip, 



where k is the wave vector of the applied plane waves and ip 
determines the wave directions. Substitute S\ and S% into the 
expressions of simulated gauge potentials A and <f> in Eq. ( TTOb 
and note that only Si and S2 are position dependent, we get 



A y = and 

A x 
hk 



- — sin 2 6cos(2(f))sinip<Jo + cos9sin(2(j))sinipa x 



A, 
hk 



— — (2 — sin 2 8) cos (2(f)) simper z , 

— (2 — sin 2 9)cosip(To + — sin 2 8cosipa z , 



(12) 



where <r x , a y , a z are three Pauli matrices and <7q is the 2x2 
identity matrix. If we set tj> = tt/4 and substitute A above 
into H e j f then we can get a spin-orbit coupling proportional 
to cos9sinipp x a x + ^sin 2 9cospp z a z which will be propor- 
tional p x a x + p z a z when we choose proper 9 and p and can 
be further turned into Rashba or Dresselhaus spin-orbit cou- 
pling through changing spin basis. However it will be a dif- 
ferent story when the Zeeman term in H e f f due to the non- 
degeneracy of the two lowest energy dressed states is consid- 
ered. 

For simplicity we assume the external trapping potential 
V = 0, so that the Zeeman term V in the simulated Hamilto- 
nian H e ff in Eq. ([8]) will be 



V 



-HA- 



sin 



2 cos(26») 



(T + HA 



sin 



2 cos(26») 



(13) 



Recall the effective Hamiltonian H e f f is simulated under the 



conditions 



{hk) 2 



<C hA and ®_ <g; ft A. Now we show under 



2m 2m 

these conditions some terms in H e f f can be further neglected. 
First we find the potential $ satisfies 



I*, 



(hk) 2 n n sin 
< ±—^sin 2 9 < ftA 



2m 



cos(29)' 



(14) 



therefore $ can be ignored compared to V. Second we write 
cos8 in A x of Eq. (fT2b as 1 — 2sin 2 (9/2) and then substi- 
tute A x and A z into H e ff, we can get some energy terms 
proportional to sin 2 9 or 2sin 2 (9/2) 12711 . which are of mag- 
nitude mai[^-, ] and can also be ignored compared to 
V. Therefore in H e f f we can safely write $ = and 

A x = hksin(p[sin(2(f>)a x — cos(2<fi)a z ], (15) 
A z = hkcostpcjQ. 

We note that this result is obtained due to the non-degeneracy 
of the dressed states \ei) and |e2), not by assuming 8 is close 
to zero. From Eq. ( Tl3T > we find that only the movement in the 
x direction of the atom is coupled to its pseudo spin, and the 
coupled movement is governed by 



H xs = tt*Z Ax) \ \ 



2m 



2m 

+ 2ap x [cos{2(j))a z — sin(2(f>)a x 



(16) 



ha. 
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where a = -^hksirup, h = fi-^- 2 cos(2B) and a constant 
term c = sin 2 w — fcA n sm L 6 L is neglected. If we use 

2m ~ 2cos{2&) ° 

\e[) = cos(j)\ei) — sin4>\e2) and \e' 2 ) = sin<j>\ei) + cos<p\e2) 
instead of \e\) and £2) to represent spin up and spin down 
respectively, then the coupled Hamiltonian will be 

H'xs = + 2ap x <r z + h[cos(24>)a z + sin(2(f>)a x }, (17) 
zm 

where 2ap x a z represents spin-orbit coupling. Thus we have 
shown the simulated spin-orbit coupling in our concerned A 
configuration can only be of a form similar to p x a y even when 
the Rabi frequencies are not much smaller than the detuning. 

The roles of h and <fi seem clear in H' xs ; one controls the 
magnitude of the "magnetic field" and the other controls its 
direction. However, as we change <f>, which is determined by 
the relative magnitude of two lasers, not only H' xs but also 
the spin basis states \e[) and |e 2 ) will be changed, and it may 
be not easy to see the underly physics. Now we assume HA 
is so bigger that a small 8 = 9 m can lead to ha z dominat- 
ing H xs , then we can study the underly physics only varying 
6 between — 8 m and m , Since 8 is small, i.e., the Rabi fre- 
quencies are much smaller than the detuning, there is \e 2 ) = 
cos(/>e _lSl |1) + sin0e~ iS2 |2) and our simulation scheme re- 
duces to the Raman process in recent experiment II 1 111 . At this 
time we get \e[) = -e~ ?;S2 |2) and \e' 2 ) = e" iSl |l), which 
are independent of <f>. Experimentally we can study the phase 
transition of H' xs due to the change of <fi and 8 as in Refs. 
ElIII. 

IV. DISCUSSION AND SUMMARY 

We have given an example that the simulated spin-orbit 
coupling in the A configuration can only be of a form sim- 



ilar to p x a y even when the Rabi frequencies are not much 
smaller than the detuning. The same conclusion can also 
be obtained when we consider the spin-orbit coupling sim- 
ulation schemes in Refs. jzTl I22I1 . So if we want to get 
a more general spin-orbit coupling in our concerned A con- 
figuration, we should find a way to eliminate the Zeeman 
term in H e ff due to the non-degeneracy of the two lowest 
energy dressed states. If we assume the trapping potential 
V = Vi\l){l\ + V 2 \2){2\ + Vs|3)(3| with V 1 = V 2 and 
V3 = V\ — E 2 j ' sin 2 8, then this Zeeman term will be elim- 
inated. However this trapping potential will give a coupling 
between vf^ and ^3 that cannot be ignored. How to effec- 
tively eliminate the Zeeman term due to the non-degeneracy 
of the two lowest energy dressed states is still under investi- 
gation. 



In summary, we have given an analytical expression of the 
simulated non-abelian gauge potentials in our concerned A 
configuration based on a special parameterization of the two 
Rabi frequencies. We have shown the simulated spin-orbit 
coupling in our concerned A configuration can only be of a 
form similar to p x a y even when the Rabi frequencies are not 
much smaller than the detuning. 
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